Electron transport properties of MgB2 in the normal state 
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Abstract 



We have measured the resistivity and the Seebeck coefficient of a MgB2 sin- 
tered sample. The temperature dependence of resistivity is fitted well by a 
generalized Bloch-Griineisen equation with a Debye temperature Qr of 1050 
K. The Seebeck coefficient is given by the sum of a diffusive and a phonon 
drag term and the behavior in the temperature region T c < T < 0.10^ follows 
a relationship AT + BT^ where the two terms are proportional to the electron 
and to the phonon specific heat, respectively. The phonon drag term, here 
emphasized for the first time, is rather large, indicating a strong electron- 
phonon interaction. The diffusive term is positive and increases with Al dop- 
ing. The comparison of the experimental values with calculations including 
precise electronic structure suggests that a bands give the main contribution 
to the Seebeck effect. 

The discovery of 40 K superconductivity in MgB2 [0 has stimulated a large discussion 
on the nature of the pairing and many evidences suggest a BCS-type mechanism: the isotope 
effect on T c [|], energy gap and specific heat measurements || and a negative pressure 
coefficient of T c f§ @. 

Electron transport properties may give insight into the normal state conduction process 
and on the electronic structure, to understand whether MgB2 can be considered a "simple" 
metal in which electron-phonon interactions play the most important role, or electron cor- 
relation has to be taken into account. From this point of view some measurements indicate 
that MgB 2 behaves indeed like a metal: magnetoresistivity is large and follows a generalized 
Kohler's rule @ , the Seebeck effect || f| |§ is small, positive and nearly linear; on the other 
hand, some Hall effect measurements || look very similar to those of cuprates decreasing as 
1/T with temperature. Therefore, for a better understanding of basic properties of MgB 2 , 
new measurements are required possibly on high quality samples because, in sintered sam- 
ples, the presence of grain boundaries can affect the transport properties. This is not the 
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case of the Seebeck effect, which, as far as the scattering processes at the grain boundaries 
may be considered elastic, is not affected by granularity, as well proved in cuprates [fT0|| ; thus, 
measurements of thermoelectric power (TEP) S can be an useful tool to provide information 
on the electronic structure of MgB2, as shown in ref. |J and ||, where a dependence of S 
on the pressure and on Al doping was observed; these authors, however, did not consider 
the phonon drag contribution, which in the case of strong electron-phonon coupling has to 
be important. 

In this paper we present resistivity and Seebeck effect measurements on a sintered MgB2 
sample. The compound MgB2 was prepared by direct synthesis from the elements: Mg 
and crystalline B were well mixed together and closed by arc welding under pure argon 
into outgassed Ta crucibles which were then sealed in quartz ampoules under vacuum. The 
samples were slowly heated up to 950 °C and maintained at this temperature for 1 day. 
X-ray diffraction shows values of the lattice parameters, a = 3.087(1), c = 3.526(1) A, and 
the absence of extra reflections. The specimen for transport measurements, (2 x 3 x 12 
mm 3 ), has been prepared by pressing the powders in a stainless-steel die into a pellet which 
was then sintered by heat treatment at 1000 °C for 2 days. 

The resistivity measurements were performed using a standard four probe technique and 
the TEP was measured using an a.c. technique described elsewhere |ll| with sensitivity of 



0.5% and accuracy of 1.5%. The gradient applied to the sample was varied from 1 to 3 
K/cm and the frequency from 0.003 to 0.008 Hz; the data were acquired with a slowly 
rising temperature (1 mK/ sec). 

The resistivity measurements are presented in Fig. 1 up to 300 K; the transition region 
is enlarged in the inset. The critical temperature defined at half of the transition is T c = 38 
K with amplitude AT C ~ 0.3 K. The room temperature resistivity has a value of 130 pflcm, 
whereas the resistivity at T = 40 K is 40 fiQcm. 

We fit the temperature dependence of the normal state resistivity to the expression: 

p(T)= Po + p ph (T) } (1) 
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where po is the temperature- independent residual resistivity and p p h(T) the phonon- 
scattering contribution assumed of the generalized Bloch-Griineisen form: 

Pph( T)^ m -^B„{ljJ Jr -f^— y (2, 

where Qr is the Debye temperature, p' is the temperature coefficient of resistivity for 
T 3> Or and m = 3—5. Eq. (2) reduces to p(T) = p +constxT m for T <C Or] indeed we see 
that from 40 to 100 K the resistivity is well fitted by a power law p(T) = po + const x T 3 . 
The best fit to our data is obtained with m = 3, p — 39.7 pflcm, Qr = 1050 K, p' = 0.49 
pflcm/K and is shown in Fig.l as a continuous line. The Qr value is in fair agreement 



with 0£> obtained from heat capacity measurements although lower values have also 
been reported M. p' which is proportional to the dimensionless Xt r coupling coefficient, has 
the same value that in A15 |fT3|, consistent with a moderately strong coupling. The low 



temperature T 3 behavior, common in transition metals where is related to inter-band 
scattering processes, suggests that also in MgB 2 the main contribution to the electrical 
conductivity comes from a more mobile band, whose carriers are scattered into empty states 
of another less mobile band. 

The TEP measurements are shown in Fig. 2, where we see the transition around 37 
K and a continuous increase from 45 K, with curvature changing from positive to negative 
above 150 K. Above the transition region the TEP appears in excellent agreement with the 
literature data [|] || both in value and in behavior. This contrasts with the wide spread 
of resistivity data available in the literature, which indicates a smaller influence of disorder 
or granularity on TEP. 

The observed temperature behavior can be analyzed considering first the diffusive 

contribution to S given by the Mott formula: 

*=£^, (3) 
Sea 

where e is charge of carriers, a is the electrical conductivity and a' = J^crfe) \ £F . Note 
that Ep must be counted upward for electrons and downward for holes. In the isotropic case 



and if the relaxation time r is independent of the energy (this is the case for scattering with 
grain boundaries), a' J a = 3/2ep, independent of the scattering processes and equation (3) 
becomes: 



n 2 K 2 T 



2e e F ' V ; 

Second, we must consider a phonon drag term S g , arising for temperatures lower than the 
Debye temperature, when the phonon relaxation time for interaction with other phonons 
and impurities is much longer than the relaxation time for phonon-electron interactions. 
For pure isotropic metals and considering only electron-phonon normal processes, an upper 



bound for S g can be estimated to be of the order [15 



where C p h is the phonon specific heat per unit volume. Thus, for T -C 0d , and assuming 
S g to really reach this upper bound, the total TEP 5* will take the form: 

S = S d + S g = AT + BT 3 , (6) 

where: 

A = ^— 7) 

2e e F K J 



B = A = £5±^J_. (8) 

3ne e n a 5 9^ 



Here (3^ = 9N J[ B is the coefficient of the low temperature phonon specific heat, iV is 
the number of atoms per unit volume, and n a is the number of valence electrons which is 
4/3 in MgB 2 . Equation (6) can be compared with the experimental TEP of MgB 2 in the 
temperature range Tq < T < CIO/) ~ 1007^ which overlaps for a large amount with our 
measurements range. 

The inset of Fig. 2 shows the ratio S/T as a function of T 2 for 40 K < T < 100 K. The 
data show a linear behavior up to a T 2 value of 8000 (T = 90 K) and then begin to bend. 



The best fit performed in the range 2000 K 2 < T 2 < 8000 K 2 is plotted as a continuous line 
and the fit parameters are A = 1.76 x 10" 2 fiV/K 2 and B = 1.26 x 10~ 6 fiV/K 4 . 

In Fig. 2 eq. (6) with A and B given by the fit is plotted as a continuous line, while 
the diffusive term Sd = AT is plotted as a dashed line. The experimental curve is well 
fitted by eq. (6) up to 100 K, above which the data change curvature and tend to increase 
linearly, with nearly the same slope of Sd- The experimental phonon drag term defined as 
S — Sd, plotted as a dotted line, tends to saturate above 250 K. This is exactly what it is 
expected for S g ; in fact, increasing the temperature, phonon- impurity and phonon-phonon 
processes become more important, and the phonon drag falls, causing a peak in the TEP. 
In metals as Cu, Ag, Au, Al, the phonon drag peak occurs at about 0_d/5 [p~5|| , while in our 
case the peak is not yet reached at 0£>/4 giving a further evidence of the importance of the 
electron-phonon coupling in MgB 2 . 

We therefore find that the diffusive and phonon drag terms contribute nearly equally to 
the TEP. This was not recognized both in ref. || and || because only a diffusive term, even 
not vanishing at T = 0, was considered in order to fit the data. 

To further verify the consistency of the model, let us try now to relate the coefficient A 
and B with some microscopic parameters. We first start neglecting the multi-band character 
of MgB2 and we pursue in our naive model of isotropic free electrons. Its reliability will be 
discussed afterwards in the light of the band structure effects. 

The coefficients A and B of our sample, and those extrapolated from data in literature 
for pure || and Al doped MgB 2 |§, are summarized in table I. The Fermi energy Ep and the 
Debye temperature Qp> are obtained directly from A and B (eq.s (7) and (8)): the Debye 
temperature of about 1400 K, is only the 30% higher than 0# and Ef for pure MgB2 is of 
the order of 2 eV and becomes 1.2 eV for Al doped samples: since the Al doping raises the 
Fermi level, the decreasing of Ep with Al doping is a further evidence (in addition to the 
positive sign) that Sd is dominated by holes. 

Let us now turn to the TEP as resulting from the band structure of MgB 2 . Two types 



bands contribute to the conduction |L7[ [|lql : two a bands, deriving from the p xy states of 



B and two tc bands deriving from the p z states. The two sets of bands have very different 
dimensional character, the a bands being of hole-type and nearly 2D, and the 3D tt bands 
mainly of electron-type. Great relevance in the discussion of the pairing mechanism has 



been given to the a bands |Tj|, and the positive sign of the Hall coefficient and of the 
TEP, as well as the increasing of the latter with Al doping, confirm the importance of these 
bands in the transport properties. Starting from the precise electronic structure calculations 



described in Ref . |L7| , we have computed the Seebeck tensor Sd , as a function of the chemical 
potential shift, to reproduce within a rigid band scheme the Al doping. These calculations 
are performed using the the scheme described, e.g., in ref. |17|,[19|. Because of hexagonal 
symmetry, the independent components are S xx = S yy and S zz . Briefly, if Q is the unit cell 
volume and fj(k, n) are the cartesian components of the Fermi velocities for the n — th band, 



we obtain the TEP components Sa as: 



e fdeau(e)(-- d£/ 



SdT) = ^ J Z*°.,Z ( 9 ) 



where 



e 2 

aij(e) = — ^Vi(k,n)vj(k,n)T(k,n)8(e(k,n) - e) (10) 

k,n 

If we consider r to be isotropic and independent of energy, it will not affect the final result. 



As pointed out in Ref. [19], Su vanishes in the approximation —j- — 5(e — /i), and it would 



therefore be appropriate to include the energy dependence of r as well. Unfortunately, 
it is not easy to obtain a consistent r(e) approximation, and we therefore use the most 
conservative, r = const, approach. 

Sd shows the expected linear behavior as a function of temperature, and we therefore plot 
in Fig. 3 Sd/T , as a function of the electron doping, in a rigid band scheme. Fig. 3(a) gives 
the tensor components, while to compare with experiment we show in Fig. 3(b) the average 
of Sd/T over directions, Sd/T . Since both numerator and denominator in the definition 
of Sd/T contain a band summation, there is no clear cut distinction between the a and ir 
contributions (in fact, when we consider the a bands alone, the resulting Sd/T reproduces 



the free-electron results along the x and y directions, and is zero along z, up to the 2D to 3D 
crossover when \x crosses the a band top at T, and S zz grows to finally give isotropic results). 
With this warning, in order to better understand our results, we decompose Sd/T in terms of 
a and tt bands contributions (i.e. we decompose the numerator of eq. (9)). We notice that 
for small values of doping the dominant, positive contribution comes from a bands; when the 
chemical potential goes beyond the A point maximum this contribution disappears and the 
resulting S&/T is much smaller and negative, identical to the 7r contribution. The smallness 
of the latter relative to the a contribution can be easily understood: the numerator of eq. 
(9) essentially monitors the on{e) derivative with respect to e, and this quantity is much 
larger for the a bands. As a function of the doping x, Sd/T initially grows, and then it 
bends down, especially when \x goes beyond the 2D to 3D crossover, but also because of the 
change the 7r contribution. 

The calculated values are smaller, by a factor of about 1.7-2, than the experimental 
values of table I. The trend as a function of doping, however, is similar in experiment 
and theory; it would be interesting, in this respect, to obtain the experimental values for 
larger doping. While the agreement between theory and experiment is not quantitatively 
good, it may considered to be satisfactory based on the fact that our calculations do not 
include the energy dependence of r, which may result into quantitative differences As 



a speculation, we may say the following: our understanding of MgB 2 indicates a strong 
electron-phonon coupling for the a band top (mostly with the E 2g phonon mode). Very 
likely, such a strong coupling will depend on the energy location, relative to the a band 
top, resulting into an energy dependent r(e). We also stress that our calculations do not 
contain any form of renormalization; it has been a controversial question, in the literature 
|2T)|| , whether renormalization should affect the bare band-strutture results. If this would be 
the case, the strong electron-phonon coupling in MgB 2 would bring the theoretical value into 
agreement with experiment. We may speculate that the presence of contributions to TEP 
from different bands, having different electron-phonon couplings, should prevent a complete 
cancellation of renormalization effects in the final result. 
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In summary, we analyzed in detail the resistivity and the thermoelectric power in the 
normal state and we showed that these transport properties can described within an inde- 
pendent electron framework, by taking into account the high phonon frequencies and the 
strong electron-phonon coupling. In particular the TEP is the sum of a diffusive and a 
phonon drag terms which contribute nearly equally to it. The phonon drag term was not 
previously recognized, in fact its peak is shifted above room temperature by the high Debye 
temperature and the strong electron-phonon interaction. The diffusive term is positive and 
increases with Al doping. The comparison of the experimental values with the theoretical 
ones, calculated starting from the precise electronic structure suggests that a bands give the 
main contribution to the TEP. Further investigation will be necessary to verify this result; 
in particular transport properties on samples with higher level of Al doping will be a useful 
tool to better investigate the role of the a bands, whose relevance in the pairing mechanism 
has been strongly advocated. 
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Figure captions 

Figure 1: p as a function of temperature; the transition region is enlarged in the inset. 
The best fitting curve obtained with p = 39.7 pQcm, Qr = 1050 K, p' = 4.95 x 10 _1 
pQcm/K is reported as a continuous line. 

Figure 2: S as a function of temperature; the continuous line is given by eq. (6) with 
A = 1.76 x 10" 2 pV/K 2 and B = 1.26 x 10~ 6 pV/K 4 ; the dashed line is the diffusive term 
Sd = AT; the dotted line is the experimental phonon drag term defined as S — Sd- In the 
inset it is shown S/T as a function of T 2 ; the continuous line is the best fit performed in 
the range 2000 K 2 < T 2 < 8000 K 2 . 

Figure 3: S^/T as a function of doping, estimated in a rigid band scheme from the 
integrated DOS. (a) Tensor components and average; (b) Average Sd/T and its band de- 
composition, as explained in the text. 

Table captions 



11 



Table 1: The coefficienta A and B obtained by fitting S/T vs. T 2 (eq.(6)); Go as 
obtained by eq. (8); ep as obtained by eq. (7). 



Table I 

A (j^V/K 2 ) B (/jV/K 4 ) Q D ( K) e F (eV) 



MgB 2 


1.76x10' 


-2 


1.26x10' 


-6 


1430 


2.1 


( a )MgB 2 


1.70x10" 


-2 


1.30x10" 


-6 


1430 


2.2 


( fo )MgB 2 


2.00x10- 


-2 


1.30x10" 


-6 


1430 


1.8 


( b )Mg .95Al .o 5 B 2 


2.88x10- 


-2 


1.41x10- 


-6 


1390 


1.3 


WMgo.gAlo.iB, 


2.99x10- 


-2 


1.49x10" 


-6 


1360 


1.2 



Wref. 0; ( 6 )ref. [|. 
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